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ABSTRACT 

Context. Both superbursters and soft X-ray transients probe the process of deep crustal heating in compact stars. It was recently shown that the 
transfer of matter from crust to core in a strange star can heat the crust and ignite superbursts provided certain constraints on the strange quark 
matter equation of state are fulfilled. 

■ Aims. We derive corresponding constraints on the equation of state for soft X-ray transients assuming their quiescent emission is powered in 
' the same way, and further discuss the time dependence of this heating mechanism in transient systems. 

Methods. We approach this using a simple parametrized model for deep crustal heating in strange stars assuming slow neutrino cooling in the 
core and blackbody photon emission from the surface. 

Results. The constraints derived for hot frequently accreting soft X-ray transients are always consistent with those for superbursters. The 

■ colder sources are consistent for low values of the quark matter binding energy, heat conductivity and neutrino emissivity. The heating 
' mechanism is very time dependent which may help explain cold sources with long recurrence times. Thus deep crustal heating in strange stars 

can provide a consistent explanation for superbursters and soft X-ray transients. 

Key words. Stars: neutron ~ Dense matter - Equation of state - X-rays: binaries - X-rays: bursts - X-rays: individuals: Aql X-1, Cen X-4, 
SAXJ1808.4-3658 

1. Introduction 



Strange stars are a class of compact stars made entirely of 
absolutely stable strange quark matter except for a possible 
thin crust of ordinary nuclear matter below the neutron drip 
density set^arated from the core by a strong Coulomb barrier 
llltohl ll97rt iBav m & ChiiJll976t IWittenlll984l: iHaensel et alJ 
ll986H Alcock et alJ ll986l) . The stabihty of strange quark mat- 
ter depends on poorly constrained strong interaction proper- 
ties and remains to be decided by observ ation o r experiment 
(see Madsen ( 1999); Weber (2005); Weber et alJ (E006) for re- 
views). If true, the strange matter hypothesis would imply that 
all compact stars are likely to be strange stars (because the 
galactic flux of quark matter lumps from strange star collisions 
will trigger the hadron-to-quark transition in ne utron stars and 
supernova cores llMadsen..l98& iCaldwell & Friedman. ! 99 1.) ). 
and so strange star models should be able to reproduce ob- 
served characteristics of compact stars. Specifically strange 
stars should be able to reproduce the bursting behavior of low 
mass X-ray binaries (LMXBs) expected to contain some of 
the heaviest compact stars due to accretion, and here we test 
the consistency of the recent model for s uperburst ignition in 
strange stars bv lPage & Cummin d (l2005h with similar consid- 



erations for the quiescent emission from soft X-ray transients; 
two of the most important LMXB activities. 

Superbursts, first discovered by ICornehsse et alJ ll200ril) . 
are a rare kind of type I X-ray bursts from LMXBs dis- 
tinguishing themselves by extreme energies (~ 10"*^ erg), 
long durati ons (4-14 hours) and long recuiTence times (~ 1 
year) - see iKuulkerj ( l2004l) for a recent review of their ob- 
servational properties. The mechanism behind superbursts is 
thought to be unstable thermonuclear burning of carbon in 
the crust at densities around 1 0^ g cm~^ fcumming & Bildstetil 
2001; Strohmaver & Brown 2 0o3lCumm ing & Macbeth 2004; 
Brown 2004; in't Zand e t alJl2004 ICooper & Naravanll200l[ 
.Cumming et aL.200 5). The fuel for such an event is provided 
by the burning of accreted hydrogen and helium at the surface, 
which may take place either stably or unstably. In the latter case 
it gives rise to ordinary type I X-ray bursts. The ashes from 
this process, containing about 10% carbon by mass, are com- 
pressed under the increasing weight of the layers above, and 
when enough carbon has been accumulated it may be ignited 
and burn unstably. Observed superburst light curves require the 
ignition to take place at temperatures around 6 x 10^ K and 
column depths of approximately 10'^ g cm"^, and as the heat 
released by accretion and burning at the surface is quickly rera- 
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dialed and does not penetrate to this depth, these ignition con- 
ditions must be set by the process of deep crustal heati ng and 
therefore probe the physics of the deep crust and core (iBrowrJ 
[2004; Cooper & Naravan 2005). 

In neutron stars deep crustal heating refers to the heat re- 
leased as the sinking ashes undergo electron captures and most 
importantly - at densities above lO'^ g c m '' - py cnonuclear 

j| ll99dl2003HBrc 



reactions JHaensel & Zduniklll99(l l2003t iBrown et al.Hl998l) . 

This amounts to about 1 .4 MeV per nucleon and heats the core 
to ~ 10*^ K depending on the accretion rate. It has however been 
difficult to achieve superburst ignition at the appropriate col- 
umn depth, and in a recent investigation^C ummin^ et al](l2005l) 
showed that neutrino emission from Cooper pairing of neutrons 
at densities in the range lO'^-lO'^g cm"-' limits the crust tem- 
perature below 5 X 10^ K. Thus unstable carbon ignition seems 
impossible at the observed depth unless Cooper pair emission 
is less efficient than currently thought or some additional heat- 
ing source can be found. 

In strange stars there can be no pycnonuclear reactions be- 
cause the crust density is limited by the neutron drip density, 
but for the same reason there can be no Cooper pair emission 
from neutrons either Furthermore the energy released by mat- 
ter tunneling through the Coulomb barrier at the base of the 
crust and converting to quark matter could be much greater 
than that released by pycnonuclear reactions in neutron stars. 
Thus it seems reasonable to expect that strange stars in LMXBs 
can be heated to higher temperatures than neutron stars, and 
[Page & Cumming (2005) recently showed that this can lead 
to superburst ignition at the observed column depths if certain 
constraints on the strange quark matter core parameters are ful- 
filled. 

Soft X-ray transients are LMXBs which undergo accretion 
outbursts lasting days to years separated by long periods of qui- 
escence in which no or very little accretion t akes place and 

1032. 



the X - ray luminosity drops to ~ 10^'^ erg s ' dCampana et al.l 
Il998t iRutledge etal]ll999t iTomsick et al']l2004 . When ac- 
cretion stops, the heat deposited at the surface quickly radi- 
ates away, and so in the absence of any residual accretion 
the surface temperature should reflect the temperature of the 
core including the effects of deep crustal heating. This has 
been thoroughly investigated for neutron s tars in th e literature 
jBrown et alJl998|IUshomirsky & Ru tledge 2001; Col pi et alJ 
| 200l[lBrown et alJ2002HRutledge et al. 2002a; Yakovlev etaD 
l2003[ l2004 and references therein) establishing deep crustal 
heating as an unavoidable heating mechanism in the core, 
which should at least provide a "rock bottom" luminosity for 
these sources. Residual accretion may still play a role however 

- indeed the power law components in the spectra of many of 
these sources are often taken as a sign of residual accretion or 
shock emission from a rotation powered pulsar - and further- 
more the properties of the core remains uncertain, so the precise 
nature of the quiescent emission is not firmly established. 

Considering the many possible phases and poorly con- 
strained properties of strange quark matter one would expect - 
for the same general reasons they were able to fit superbursters 

- that a wide range of strange star models can reproduce the 
observed properties of soft X-ray transients. Our approach is 
therefore to use a simple parametrized model for deep crustal 



heating in strange stars to discuss the consistency between con- 
straints on strange quark matter derived from superbursters and 
soft X-ray transients, as such a mutual fit is likely to constrain 
the strange quark matter equation of state further 

We describe the overall structure of strange stars and their 
crusts in Sect.|2]and then turn to deep crustal heating in Sect.|3l 
In Sect.0]we compare models with different assumptions about 
the cooling mechanism to observations of soft X-ray transients, 
and we discuss the importance of time dependence in Sect. |5] 
before concluding in Sect.|6l 

2. Strange stars and their crusts 

If strange stars are stable they contain roughly equal num- 
bers of up, down and strange quarks in their core, but due 
to the higher mass of the negative s-quark, these will nor- 
mally be slightly disfavored giving the core an overall posi- 
tive quark charge to be compensated by electrons. Even color- 
flavor locked^trange quark matter, which is strictly charge neu- 
tral in bulk feaiagopal & Wilc zek 2001 ) has an overall positive 
quark charge due to surface effects ( Madsen l2000t EoOU lUsovl 
2004). Because strange quark matter is self bound the density is 
largely constant in the core dropping by only a factor of a few 
from the center to the surface, where it goes from above nu- 

Alcock et alJ 



(1986); Haenseletal. 


( 1986); Glendenning & Weber 


l|l992 


Kettneretal. (1995); 


Weber (1999); 


Glendenningl 


ll200nl 



Zdunik etal. (2001); Zdunik (2002) for details on flie global 
properties of strange stars. Since the electrostatic force is much 
weaker than the strong force confining the quarks, some of 
the electrons necessary to create an overall charge neutral ob- 
ject will form a thin atmosphere outside the quark phase held 
in place by a huge electric field. The electric field - up to 
10"* V cm"' - is capable of sustaining a nuclear matter crust 
dec oupled from the quark phase by a small electron filled 
gap d Alcock et al.ll986tlKettner et al.ll99llSteiner & MadsenI 
2005). The crust is electrically neutral in bulk however, so be- 
yond the boundary layer its structure is given by hydrostatic 
equilibrium with an equation of state corresponding to the 
chemical composition of the ashes from the surf ace burning 
dGlendenning & Weber[| l992": 'Zdunik et^al.ll200ll) . The crust 
(Izdi 



mass roughly follows 



Munik etalj2001 ) 



Mc 



GM 



2GM\ 



8.1 X 10"^— i-^ 



1 



0.295^ 



(1) 

(2) 



where Ph is the pressure at the bottom of the crust, M, = M/Mq 
is the stellar mass in units of the solar mass, R(, = R/10 km is 
the stellar radius, pb is the density at the bottom of the crust, 
pu - 7.8 X lO" g cm"^ is the neutron drip density, y ^ 4/3 is 
the adiabatic index in the crust and the last equality assu mes the 
crust composition found by Haensel & ZdunikI j2003h for ac- 
creting neutron stars. The density at the bottom of the crust may 
be limited by the neutron drip density if the Coulomb barrier 
from the core is sufficiently high, which depends on the details 
of the equation of state of strange quark matter. In most cases 
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this will not be possible though and the maximum crust density 
is then determined by the rate at which ions tunnel through the 
Coulomb barrier such that the tunneling rate for an accreting 
star is equal to the accretion rate or corresponds to th£ inverse 
stellar age if the star is isolated (see lSteiner & MadsenI fcoOSi) 
for a detailed discussion). The tunneling rate close to this equi- 
librium is very sensitive to the crust density so there is no great 
difference between these two cases in terms of crust density, 
but for the same reason one may expect that in soft X-ray tran- 
sients the tunneling rate will vary with the accretion rate - we 
return to this point in Sect.|5] 

Since the crust has no significant heating sources of its own 
(electron captures in the criist am ount to less than 40 keV per 
nucleon faaensel & Zdunikl2003l) ') its thermal profile is set by 
the diffusion of heat from the core towards the surface. The 
surface temperature, Ts, is therefore only a function of the core 
surface temperature, Ti, and thus probes the process of deep 
crustal heating. A detailed account for such 7$ - Ti relations 
can be found in e.g. Yakovlev et al. (2 004i), but in a r eason- 
able approximation T^g ^ O.IT^^ dOudmundsson et aDfl983l) . 

which will be sufficient for our purpose. Here Ts.e = Ts/lO^ 
K and T^g = Ti/lO'' K is actually the temperature at a den- 
sity of 10'*^ g cm""*, but as the crust is largely isothermal at 
this depth, one can take this to b e the core surface temperature 
Joudmundsson et al, 1983; Yakovlev et alj.2 004). 

In the discussion above and in what follows we have as- 
sumed the conventional description of strange stars with a 
sharp core surface, an electron filled gap and a nuclear cr ust. 
It should be mentioned however that Jaik umar et alJ ( l2006l) re- 
cently showed, that if the surface tension of strange quark mat- 
ter is not too high, it may be unstable to phase separation at 
low pressure leading to a crust of quark nuggets embedded in a 
degenerate electron gas, which would in many ways resemble 
the nuclear crust on neutron stars. The consequences of such 
a crust in an accreting system are not clear at present, and we 
restrict ourselves to the traditional scenario. 



3. Deep crustal heating in accreting strange stars 

De ep crustal heating in str ange stars was recently studied 
by IPage & Cummind ( 12005 ) as a means of achieving su- 
perburst ignition conditions at the appropriate column depth 
and thus constrain the properties of strange quark matter 
(see also the work on continuously accreting strange stars by 
[Miralda-Escude et al. ( 1990)). It was shown that observed ig- 
nition conditions are achieved for a wide range of crust density, 
Pb, strange quark matter binding energy, Qsqm, neutrino emis- 
sivity, ey, and thermal conductivity, K, provided that the direct 
Urea process does not operate in the core implying that strange 
quark matter should be a color superconductor Furthermore 
a set of useful analytical approximations to the core surface 
temperature were shown to reasonably approximate numerical 
calculations. Here we will use those same approximations to 
find the corresponding constraints on strange quark matter in 
soft X-ray transients and a relationship between the observable 
parameters of superbursters and soft X-ray transients. For clar- 
ity and to emphasize the approximations we therefore briefly 



review the calculation bvlPa ge & Cummin^ ll200.'5l) in the fol- 
lowing subsection. 

3.1. Approximate thermal structure 

As accretion on the strange star drives the crust towards its 
maximum density, the rate at which nuclei (or neutrons if 
the neutron drip density can be reached) tunnel through the 
Coulomb barrier into the quark matter core will increase un- 
til an equilibrium crust density is reached. This density may 
vary as a function of the unknown Coulomb barrier height at 
the quark matter surface but this would not change the qualita- 
tive features of the model and for our purpose we will assume 
that the equilibrium crust density is ~ 10" g cm"-'. As the nu- 
clei cross into the strange quark matter core at the rate m (in 
units of nucleons cm"^ s"' unless otherwise stated) and are dis- 
solved, they release an energy, 2sqm < 100 MeV per nucleon 
(the actual value depends on uncertain values for strong inter- 
action parameters but is fixed if these parameters were known). 
This heats the core until an equilibrium or steady state temper- 
ature is reached as the heating balances energy loss by neutrino 
emission in the core. This is assumed to take place at a (slow) 
rate of 



Qv X Ta erg cm s ' , 



10 



22 



(3) 



with Tg - r/lO^K, corresponding to modified 
Urea and bremsstrahlu ng processes (Paseetal. 2005j 
lYakovlev & Pethickl2004 . We thus neglect a slow dependence 
on density in Qy which would only amount to a factor of order 
unity since the core density is nearly constant. 

We will - at first - assume that this equilibrium is simply 
determined by the average accretion rate, which would be rea- 
sonable in systems where the tunneling is more or less continu- 
ous. We return to this point in section|5]however as this steady 
state assumption may not be sufficient for all systems. 

Following Page & Cumming (2005) we model the thermal 
profile of the core in a plane parallel approximation with heat 
transport taking place by conduction and with a heat sink from 
neutrino emission. The inward directed heat flux from the con- 
version of nucleons to strange quark matter, F, is determined 
by 



dF 



-pK 



.dT 
dy 



(4) 



where p is the density, K is the thermal conductivity of strange 
quark matter, assumed to lie the range 10''^ - 10^^ cgs, and y is 
the column density. Again we will neglect slow dependencies 
of K on density and temperature so assuming that p, Qy and K 
are constant in the core this gives 



T ^a{y + by 



(-2/7) 



11 7 
a - \ — — 



1 



Qv 



2 9(10'^fp^-K 



-1/7 



(5) 



where b is a constant of integration determined by the require- 
ment that all heat entering the core is emitted as neutrinos, 
riiQsQM - P ' ^ eydy. Assuming that the integral is dominated 
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by the surface contribution this gives for the core surface tem- 
perature 

where 2^,21 = Qv/iO^\ K20 = /T/IO^" cgs and mEdd - 1-6 x 
10"^ Mo yr^^ I AnR^niu is the Eddington accretion rate in units 
of nucleons cm"^ s"' with R being the stars radius and m„ the 
atomic mass unit. 

If however the thermal conductivity of strange quark mat- 
ter is sufficiently large the core becomes isothermal with the 
temperature given by ^JiR^mQsqyi - i'^^R^ /^)QvT^g so 

where is the core radius in units of 10 km. Equating 
this to Eq. (|6j shows that the core becomes isother- 
mal at a critical thermal conductivity Ka-n - 4.3 x 
lO^i cgs e;^^2'i(GsQM/100MeV)^/^(m/0.3 mEdd)'^'- 

It should be noted that although Eqs. ^ and Q are cer- 
tainly oversimplified they were shown by Page & Gumming, 
(l2005h to agree reasonably well with numerical calculations of 
a full general relativistic model. Furthermore they offer a phys- 
ical transparency, which given the many poorly constrained 
properties of strange quark matter, is perhaps just as desirable 
in what follows as full numerical models. Specifically Eqs. (|6jl 
and are general enough not to depend on the properties of a 
particular strange quark matter phase - except for the assump- 
tion of slow neutrino emissivity necessary to obtain superburst 
ignition. They are in fact general enough to also describe neu- 
tron stars at the same level of detail, and indeed Eq. d lt corre- 
sponds entirely to case (ii.a) in Yakovlev et al. '('2003') for low- 
mass neutron stars with slow neutrino emission. Thus much of 
what follows would also apply for neutron stars - although in 
that case it has of course been done in much greater detail. 

In deriving Eqs. (|6} and we assumed that cooling by 
neutrino emission dominated photon emission from the sur- 
face. This is certainly the case in superbursters, where the 
core is very hot, but for cold enough cores this may not be 
the case and one should take the surface luminosity into ac- 
count. Approximating the surface with a blackbody we can es- 
timate when this becomes important by setting , AnR^i nQspM = 
AnR^crT^, which using the lOudmundsson et alJ ill 9831) relation. 
Ti g ^ 0. 1 g, gives a core temperature 

Equating this to Eq. shows that cooling by photon emission 
becomes important in the isothermal case below an accretion 
rate of 

with a similar expression giving a somewhat lower limit for a 
non-isothe rmal core. Simil ar conclusions for neutron stars were 
reached bv lYakovlev et alJ (EoO^ . Given the right combination 



of core parameters photon emission may therefore be important 
in soft X-ray transients with the lowest accretion rates as low 
as a few times 10"'^ Mg yr"'. 

In the following subsections we shall ignore photon cool- 
ing from the surface, which allows analytical constraints on 
the core parameters and simple scaling relations for the sur- 
face temperature independent of assumptions about the core 
parameters. This approach is justified for high 2sqm, but we 
return to discuss the importance of photon cooling in Sect.l^as 
we shall find neutrino cooling unable to account for the coldest 
soft X-ray transients. 

3.2. Constraints on the core parameters 

Eqs. (|6j and relate quark matter properties represented by 
Qv, 6sQM and K to the stellar accretion rate and core sur- 
face temperature under the assumption that the accretion rate 
is high enough for neutrino cooling to dominate. In case there 
are independent constraints on the core surface temperature 
and the accretion rate Eqs. (|6} and can also be regarded as 
necessary relations between the quark matter core parameters 
Qv, Qsnivi and K. This a pproach was taken for superbursts in 
Pa2e & Gumming assuming that the crust below the su- 

perburst ignition depth is approximately isothermal so the core 
surface temperature can be approximated by the superburst ig- 
nition temperature, T^g ^ 0.7 at m = 0.1 - 0.3 mEdd- We will 
make the same assur nption and take over the su perburst param- 
eter constraints fromlPa ge & Gumrnin 

Assuming that soft X-ray transients are also powered by 
deep crustal heating, a different independent constraint on the 
quark core parameters can be obtained by applying Eqs. (|6} 
and to observations of the accretion rate and surface temper- 
ature of soft X-ray transients using a Ts - T, relation such as the 
one in Yakovlev et al. (2004) or Gudmundsson et al. ( 1983). 

Aql X-1 is an example of a well studied soft X-ray transient 
which has (m) 10''° y r'', ^ 1 13 eV jYakovlev et alJ 
l2003t iRutledge et alJl2002 al) and therefore T[ g ^ 0.3 using a 
canonical mass-radius relationship, R - 10 km, M/R = 0.2, 
and Ts = T^/(l -2GM/R)^'^. The relations between the quark 
matter core parameters thus obtained for superbursts and for 
the soft X-ray transient Aql X-1 are compared in Fig. [2 If we 
keep Qv fixed and increase Qsqm then K must also increase to 
produce the required temperature until the critical conductiv- 
ity is reached and the star must be isothermal. Therefore only 
the regions below the isothermal curves contain possible solu- 
tions for each class of object and there can be no common so- 
lution for superbursters and Aql X- 1 between the two isother- 
mal curves shown. It is clear that the range of quark matter 
parameters permitted for the soft X-ray transient Aql X-1 is 
significantly smaller than that allowed by the superbursters; in 
particular the quark matter binding energy per baryon, Qsqm, 
necessary to fit the soft X-ray transient must be in the lower part 
of the range allowed for fitting superbursters, other parameters 
being equal. We further note that the non-isothermal relations 
from Eq. (|6j agree remarkably well in the range where common 
parameters fit superbursters and Aql X- 1 , whereas no simulta- 
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neous isothermal solutions for the two classes of objects seem 
to exist. 




le+18 



le+22 



Fig. 1. Constraints on quark core parameters obtained from 
superbursters (solid) and the soft X-ray transient Aql X-1 
(dashed) assuming slow neutrino cooling for both. Qsqm is 
given as a function of Qy for different choices of K. Thin 
lines show non-isothermal constraints from Eq. (|6jl with K - 
W^, 10^°, 10^', 10^^ cgs from the bottom up. Non-isothermal 
constraints for Aql X-1 are nearly indistinguishable from those 
for superbursters. Thick lines show isothermal constraints from 
Eq. Q. We take m - 0.3mEdd for the superbursters. 



For a fixed (albeit unknown) set of strong interaction pa- 
rameters neither the quark matter binding energy 2sqm, nor 
the quark core surface density should vary between stars, be- 
cause both are given directly by the strong interaction parame- 
ters. The neutrino luminosity is dominated by the core surface 
contribution and strange stars are furthermore characterized by 
an almost constant density in the core unless we are very close 
to the maximum mass strange star, so the slow variations with 
density of Qy and K should be unimportant. K may be inversely 
proportional to the temperature ( HaenseL1991) and can there- 
fore vary by a small factor between the coldest soft X-ray tran- 
sients with quark core temperature ~ 10^ K and superbursters 
with core temperature ~ 7 x 10*^ K. If no density thresholds be- 
tween phases or emission mechanisms are crossed when going 
from superbursters to soft X-ray transients it should therefore 
be reasonable to expect approximate agreement between the re- 
lations in Fig.n as seen for the non-isothermal solutions, but 
not order of magnitude differences as found in the isothermal 
case. 

3.3. Scaling relations for the surface temperature 

In this subsection we will use Eqs. (|6j and Q to relate the 
parameters for strange stars with deep crustal heating that can 
fit soft X-ray transient data and superbursters simultaneously, 
assuming slow neutrino cooling to be the dominant cooling 
mechanism. In particular we will derive relations between sur- 
face temperature and accretion rate for soft X-ray transients for 
strange star model parameters that fit superbursts. These rela- 



tions will then be compared to observational data for soft X-ray 
transients in the following Section. 

To the extent that Qsqm and Qy do not vary between stars 
we can utilize Eqs. (|6} and (0 to compare the other parameters 
for two arbitrary strange stars with deep crustal heating. In the 
following we shall in particular be interested in cases where 
one of the stars is a superburster (SB) and the other a soft X- 
ray transient (X). From Eq. (|6} we have in the non-isothermal 



asB = 



p9/2 
i,9 



0.3 /WEdd 



SB 





/ 

7,9/2 
i,9 

\ 






\-^^20,SB/ 


0.3mEdd 





(10) 



and from Eq. Q in the isothermal case 

<m) 



asB = 

R(>,st 



i,9 



0.3 mEdd 



SB 



9,1 



0.3 mEdd 



(11) 



In a soft X-ray transient the temperature at the quark core 
surface can be obtained from a Ts - T, relation with Ts - 
r|°/(l - 2GM/Ry^^. The average accretion rate can be esti- 
mated from the outburst luminosity, L„ ~ (AM/to)(GM/R), 
where AM is the accreted mass and is the outburst dura- 
tion, so if tq is the time spent in quiescense (m) - j^^fjp^- 
In Eqs. (I10> and M \\ we should therefore further allow for dif- 
ferent MIR ratios. If we use a detailed Ts - T, relation such as 
in Yakovlev et al. (2004) there is no point in doing this analyt- 
ically. As noted such relations can however be approximated 
reasonably well by T/g ^ 0.1 Tg^^ which allows a physically 
transparent interpretation. 

We will use Eqs. (I10> and (IIU to compare superbursters 
and soft X-ray transients so writing a^B for the left hand side 
scaled to fit the superburst data and dropping the subscript on 
the soft X-ray transient parameters we get 



rs,6 =10 



1/2 



Q'SB 



^20,sb\'^^ {m) 



K20 



0.3 mEdd 



1/9 



Q'SB =0.2 





<m) 


-r 






0.3 mEdd 




SB 



(12) 
(13) 



in the non-isothermal case and correspondingly in the isother- 
mal case 



1/2 



Ts,6 =10 



asB =0.058 



asB 



Re.SB \ <m> 



Re j 0.3 mEdd 



1/16 



m 


(ih) 








0.3 mEdd 




SB 



(14) 



(15) 



As a last possibility we note from Fig.^that solutions can be 
found with non-isothermal superbursters and isothermal soft 
X-ray transients. In this case we use A'crit in Eq. (|6j to obtain a 
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relation between (2sqm and Qy which when inserted in Eq. Q 
gives 



Ts,6 =10'^' 
asB =0.154 



0.3 niEdd 



1/9 



((^) 




-n 






0.3mEdd 


) 


SB 



(16) 
(17) 



These expressions can be used to test the constancy of the 
quark core parameters when going from superbursters to soft 
X-ray transients, with a^B incorporating the constraints im- 
posed by superbursters. We note that small variations in the 
quantities in square brackets will be strongly suppressed and 
specifically it will not matter whether we use a canonical mass 
radius relation for neutron stars to find the accretion rate or one 
more relevant to strange stars where the radius may be smaller 
by a few kilometers. This would translate into a temperature 
correction of about one percent which is far below the obser- 
vational uncertainty so we can use the accretion rates in the 
literature even though they are derived for neutron stars. There 
is also no reason to expect that (/?6,sb/^6)'^'^ will deviate much 
from unity so we drop this in the following. The (/Tsb/^)'^'^ 
term can give a factor ~ 0.9 if K depends on the temperature 
and Q-sB :S 0.2 so we keep both of these terms but note that 
uncertainties in these parameters are strongly suppressed too. 

4. Comparison with observations 

The observational sample most commonly used to test rela- 
tions for deep crustal heating such as Eqs. MA\ and M21 con- 
sists of soft X-ray transients in the field with short outburst 
durations and recurrence times, for which X-ray observations 
in quiescense have determined or limited the surface tempera- 
ture. The outburst duration should be short so the crust has not 
been significantly heated by the release of nuclear energy from 
the accreted matter and the surface temperature is therefore set 
by deep crustal heating. The recurrence time should be short 
enough that at least a few outbursts have been observed allow- 
ing a reasonable estimate of the time averaged accretion rate 
- although one should remember that the last 37 years of X- 
ray observations are not necessarily representative of the ~ 10"* 
years of accretion history needed to reach thermal equilibrium. 
Similarly since the accretion history for stars in globular clus- 
ters is uncertain one should only use systems in the field. 

Useful lists of such systems a nd their proper t ies inc lud- 
ing reference s can be found in e.g. lYakovlev et all (12003) and 
iTomsick et al.. (.2004). We will take this data, summarized in 
Tabled as our observational basis and briefly discuss each sys- 
tem below. The best examples of such systems are Aql X-1, 
Cen X-4 and 4U 1608-522 being short outburst field soft X-ray 
transie nts with detected thermal components. iTomsick et alJ 
(120.04') discuss the outburst history of these sources and we cite 
their accretion rates which are based on the review of soft X-ray 
transient activity bv lChen etalJ(ll997l) using 



<m) = 



ec2(33 yr) 



(18) 



where Lpeak is the time averaged peak outburst luminosity, 
is the number of outbursts, / = 0.4 is a factor which accounts 
for the shape of the light curve, e - 0.2 is the fraction of rest 
mass energy released during the accretion and the time average 
is taken over 33 years of observation. Aql X-1, Cen X-4 and 
4U 1608-522 have gone though 21, 2 and 16 outburst respec- 
tively between 1969 and 2003 allowing estimates of the accre- 
tion rate. Temperatures are found from spectral fits to absorbed 
neutron star atmosphere models assuming a mass of 1 .4 M© - 
this also yields the best fit radius cited. Cen X-4 and some ob- 
servations of Aql X- 1 further require a power law component 
possibly indicating residual accretion during quiescense. 

For SAX J1808.4-3658, SAX J1810.8-2609 and XTE 
J2123-058 no thermal components have been detected and so 
we only cite upper limits to the surface temperature. SAX 
J1810.8-2609 and XTE J2123-058 have each shown only one 
outburst so the accretion rate is actually an upper Umit. The ac- 
cretion rate cited for SAX J1808.4-3658 is based on the two 
outbursts in 1996 and 1998, but further outbur st have ocurred 
in 2000, 2002 and 2005 (see IWiinandsl ( |2004 for a review). 
These systems are therefore not expected to agree well with 
the theory, but at least the upper limits to their surface temper- 
atures should not lie significantly below the prediction of deep 
crustal heating if they have accreted for long enough to have 
reached thermal equilibrium. 

XTE J 1709-267 Hes 9-10 arcmin from NGC 6293 which 
has a tidal radius of 14.2 arcmin and so may be associated with 
this globular cluster although di stance measurements are incon- 
clusive on this (.lon ker et alJ200 4i). It has shown two outbursts 
since 1995 and has a varying power law component. 

As should be clear from the discussion above the data is 
rather uncertain and should be treated with some caution. In 
particular the time averaged accretion rates are hardly more 
than estimates, and the surface temperatures, when they can be 
found, are uncertain and may not only reflect deep crustal heat- 
ing. Nevertheless deep crustal heating should provide a mini- 
mum temperature for systems that have accreted long enough 
and so we should at least test for consistency with the observa- 
tions if not perfect agreement. In Fig.|2lwe plot the data from 
Tableland compare with photon cooling curves as well as 
Eqs. ( I12l i. il4\ and il6i . that assume neutrino cooling, for val- 
ues of asB corresponding to msB =0.1 and 0.3 /WEdd- 

The overall agreement with our predictions from neutrino 
emission in superbursts is not impressive. Except Aql X-1 and 
4U 1608-522 the observations fall significantly below the pre- 
dicted temperature for neutrino cooling, and actually pure pho- 
ton cooling does better if Qsqm ^ 1 MeV although it should 
not be dominant for the high accretion rates in Aql X-1 and 4U 
1608-522. This raises two issues. 

First it is noteworthy that the two sources which fall close 
to the superburst prediction without photon cooling, are the 
hottest and most frequently accreting sources thus resembling 
superbursters the most. In turn this raises the question whether 
the accretion history may lead to heating with a pronounced 
time dependence and whether the time averaged accretion rate 
is then still a good predictor for the temperature - we discuss 
this possibility in the next section. 
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Source 


kT^ [eV] 


R 


{m) [M^yr-'] 


Reference 


Aql X-1 


nyl 


13.2 


10-'" 


Rutledse et al. (2002a); Yakovlev et al. f2003) 


Cen X-4 


76 ±7, 85±3 


12.9 


1.9 X 10-" 


Rutledee et al. (2001); Camoana et al. (2004); Tomsick et al. (2004) 


4U 1608-522 


170±30 


9.4 


6.9 X 10-" 


Rutledse et al. (1999); Tomsick et al. (2004) 


SAXJ1808.4-3658 


<54 




5 X 10-'2 


Camoana et al. (2002); Yakovlev & Petliick (2004) 


SAX Jl 8 10.8-2609 


< 72 




1.7 X 10-'^ 


Jonker et al. (2004); Tomsick et al. (2004) 


XTEJ2 123-05 8 


<66 




2.9 X 10"'^ 


Tomsick et al. (2004") 


XTE J 1709-267 


116±2 




1.2 X 10^'° 


Jonker et al. (2003. 2004) 



Table 1. Properties of the observational sample commonly used to test deep crustal heating. The temperature quoted for Aql X-1 
is the red shifted rg°. We use = T ^l(\- IGM/RYJ^ wi th the canonical M/R = 0.2. For Cen X-4 we quote the temperatures 
from both lRutledge et alJ kOOlh and lCampana et alJ ll20Q4l) . 




Accretion rate [M^^^ yr ^] Accretion rate [M^^^ yr ^] 

Fig. 2. Left: comparing cooling by a non-isothermal (dashed and dot-dashed; cf. Eq. il2\ ) and isothermal (dashed-doubledot and 
solid; cf. Eq. (II 4» core to the observations in Tabled For Cen X-4 the quoted temperature intervals overlap. Right: the "mixed" 
case with non-isothermal superbursts and isothermal soft X-ray transients (dashed and solid; cf. Eq. il6\ ) and photon cooling 
(dotted) assuming a pure blackbody with 2sqm as labelled for each curve. 



Second we see that if Qsqm is sufficiently low, photon 
emission from the surface may be very important and so we 
should combine the two cooling mech anisms. Thus we now 
demand - similarly to the procedure in lYakovlev et all ( l2003l) 
- Ldh - 47r/?^m2sQM = Ly + Ly where Ldh is the luminos- 
ity required to carry away the energy released by deep crustal 
heating and Ly = AttR^ctT^ is the usual blackbody luminos- 
ity. Following the calculation that leads to Eq. (|6} the non- 
isothermal neutrino luminosity is 

Ly = 47r7?2 . 0.3 mEdd ■ 100 MeV ■ QlaXio^ ■ (19) 

Neutrino cooling should dominate in superbursters, Ldh = Ly 
with T^^ - 0.7, so using the constraints this imposes on the 
core parameters Eq. jl9t gives 

Ly^AnRH-^\ msBesQM, (20) 

and using Ti g = 0.1 T^g the condition Ldh = 4;r^^m2sQM = 
Ly + Ly then becomes 

m 2sQM = "iSB 2SQM 7^3,6 + ^^S • ^^^^ 



Correspondingly in the isothermal case 

m Qsqm = 1 -8 x 10-'otsb QsquT^^ + (tTI . (22) 

This can be solved numerically and in Fig. |3 we show rep- 
resentative solutions for a range in 2sqm and compare with 
pure photon and neutrino cooling. We see that photon cooling 
dominates at low accretion rates and low Qsqm whereas for 
high accretion rates and Qsqm we recover the neutrino cool- 
ing curves relevant to superbursters. In particular we note that 
6sQM =0.4 MeV while still fulfilling the superburst require- 
ments at high accretion rates provides a close fit to some of 
the sources at low accretion rates that were diflicult to explain 
with neutrino cooling alone. In this case other processes must 
be heating the surfaces of Aql X-1 and 4U 1608-522 though. 

iBrown et al ] iI199R|) similarly found that deep crustal heat- 
ing of neutron stars could explain the quiescent emission from 
Cen X-4 if only 0. 1 MeV of the energy released by pycnonu- 
clear reactions durin g an outburst was cond t icted into the core 
and deposited there. lYakovlev et"aD ( l2003l '2004) instead in- 
terpreted Cen X-4 and SAX J1808.4-3658 as massive neutron 
stars with enhanced neutrino emission. In the context of strange 
stars the energy is deposited directly in the core, but no en- 
hanced neutrino cooling is required if the quark matter binding 
energy Qsqm is small. 
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Accretion rate [M^^^ yr ^] Accretion rate [M^^^ yr ^] 

Fig. 3. Left: comparing combined photon and neutrino cooling of an isothermal core (solid; cf. Eq. M2V ) to the observations in 
Table^with 2sqm as labelled for each curve . For reference we also show curves with neutrino cooling only (dashed; cf. Eq. MAM ) 
and photon cooling only (dotted). The curves assuming Qsqm = 100 MeV and neutrino cooling only are nearly identical. For 
Cen X-4 the quoted temperature intervals overlap. Right: corresponding curves for a non-isothermal core. 



5. Time dependence of the mass transfer rate 

The quiescent luminosity fro m some o f the sources dis- 
cusse d here (i.e. Aql X-1 CRutled ge et al.. 2002a). a n d Cen 
X-4 llvanParadiisetalJ Il987t ICamnana et alJ Il997l l2004 
iRutledge et al."'200l')) have been found to vary by factors of a 
few on timescales of days to years raising the question of how 
time dependent the heating mechanism is following an out- 
burst. Furthermore the time dependence of the heating mecha- 
nism may be important in sources with long recurrence times 
and help explain their relatively low temperatures reducing the 
need to keep Qsqm small. 

To fully demonstrate either of these points would require a 
detailed time dependent model for cooling strange stars, which 
is beyond our scope here. Thus the calculations below are con- 
cerned only with the time dependence of the tunneling rate - 
and hence the heating - following an outburst, but we note that 
with respect to its consequences for the surface temperature we 
may be guided by similar considerations for neutron stars. 

In neutron stars heat is released during and after an ac- 
cretion event at the depth and on the time scale (months for 
the pycnonuclear reactions) corresponding to the individual 
electron captures, neutron emissions and pycnonuclear reac- 
tions. It then takes about 1 yr for the heat from the pycnonu- 
clear reactions to diffus e to t he surface jBrown et al.lll998l 
lUsh omirskv & Rutledg3 l200ll) and since it spreads out along 
the way, it is seen there as a heat wave lasting approximately the 
therma l diffusion time to the surface. lUshomirskv & Rutledgd 
investigated the time dependence of this process for 1- 
day long outbursts, and found that energy deposited at depths, 
where the diffusion time to the surface was shorter than the out- 
burst recurrence time, could produce significant variations in 
the surface luminosity following a pattern in time correspond- 
ing to the diffusion time for each nuclear reaction. If for in- 
stance the recurrence time was 1 year and assuming a low ther- 
mal conductivity in the crust, the luminosity variations were 
a few percent of the average, and only one heat wave could 
be seen because the diffusion time from the depth of the pyc- 



nonuclear reactions was on the order of the recurrence time. If 
the recurrence time was 30 years the variations in the luminos- 
ity could be up to 30 percent for cores with standard cooling 
and a factor of 3-4 for cores with accelerated cooling (e.g. pion 
condensates). In this case two separate heat waves were seen - 
one corresponding to electron captures in the outer crust visi- 
ble a few days after the outburst and the other to the pycnonu- 
clear reactions deeper in the crust emerging about a year later. 
For higher thermal conductivity variations were larger and hap- 
pened on shorter time scales. 

By analogy if the tunneling rate at the bottom of the crust - 
and hence the heating - is sufficiently time dependent follow- 
ing an acc retion event, we would expect t he qualitative con- 
clusions in lUshomirskv & Rutledgd ( l200ll) to apply, possibly 
leading to strong variations in the surface temperature follow- 
ing an outburst as the resulting heat waves reach the surface. 
One difference may be though that in strange stars the heat is 
released at a lower density in the crust with a shorter diffusion 
time to the surface so the two heat waves may be closer in time. 

The transfer of mass from crust to core by tunneling take s 
place at a rate of dAlcock et al.ll98dlSteiner & Madsenl2'005l) 

dMc 

= -WNA^ion X V X T(pb) 

= -AnR^dpb X V X T(pb) (23) 
= -23.8^r(pb) Mo S-' 

where is the mass of ions at the bottom of the crust, A'ion 
is the number of such ions whose motion allow them to strike 
the Coulomb barrier taken to be the number within one lattice 
distance, d ~ 200 fm, pb is the density at the bottom of the 
crust, V < 1 MeV is t he fr equency of this motion and T(p) is the 
transmission rate. In ISteiner & Madseni (l2005l) we found that 
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the transmission rate can be approximated by 



Ze,p,(Rs) 

2\l/2 



1/2 



(24) 



(25) 



where Z is the ion charge number, A is the mass number, 
e<pe(Rs) is the height of the Coulomb barrier (i.e. the electric 
potential at the core surface, Rs) and A0g - (po{R) - (f>o{Rs) 
is the difference in the (Newtonian) gravitational potential be- 
tween the surface of the star and the core. We plot this be- 
haviour in Fig. m assuming ejpgjRs ) - 2 MeV and the crust 
composition in Haens el & Zdunitd (l2003l) ("'^Ge at maximum 
density). From Eq. j23> a star accreting 10"'° Mq yr"' must 
have log t < -42 in equilibrium, so as previously mentioned 
the tunneling rate close to equilibrium is very sensitive to the 
density. 




le-Oe le-05 0.0001 0.001 0.01 

Fig, 4, Plot of Eq. (|24j assuming e^Rs) = 20 MeV and '°''Ge 
at the crust boundary. Note that for this barrier height the gap 
width goes to zero and r to 1 somewhat below the neutron drip 
density, which can then not be reached. 

Eq. Q gives a relation between Mc and pb so we can solve 
Eq. ( I23> to find the crust mass and tunneling rate for given ac- 
cretion scenarios. Assuming the same parameters as in Fig. ^ 
we show a solution in Fig. [S] with 10 day long accretion out- 
bursts separated by quiescent intervals of 30 years and an av- 
erage accretion rate of 1.9 x 10"" Mq yr"' - roughly corre- 
sponding to Cen X-4. We start the integration with a crust 
mass sufficiently low that the tunneling rate is negligible and 
let it build up until the mass transferred to the core during 
each cycle matches the average accretion rate. At this point the 
crust mass is at its maximum - as determined by the choice of 
e(pe(Rs) - 20 MeV - during the accretion outbursts and then 
relaxes to a state with very little tunneling within about a year 
after each outburst. Fig. |6l shows the corresponding result for 
a system resembling Aql X-1 accreting 10 in every 220 days 



with an average accretion rate of 10"'" Mq yr"' . Here we more 
clearly see the crust reach a maximum mass during the 10 day 
long outbursts. 



1.2997 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



„ 1.2996 
'o 

I L2995 
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3630 36-IO 3650 3660 3670 3680 36^0 

I I 



Time [yr] 

Fig. 5. A crude model of Cen X-4. Crust mass in a system ac- 
creting with 10 day long accretion outbursts separated by qui- 
escent intervals of 30 years and an average accretion rate of 
1.9 X 10"" Moyr"'. The inset shows a few equilibrium cycles. 
The peaks seen here have a flat top which is not resolved in 
time on the scale shown. 




Time [yr] 

Fig. 6. A crude model of Aql X-1. Crust mass in a system ac- 
creting with 10 day long accretion outbursts separated by qui- 
escent intervals of 220 days and an average accretion rate of 
10"'*' Moyr"' . The inset shows a few equilibrium cycles. 

Thus practically all the heating takes place during and im- 
mediately after the outburst when the transfer rate is essentially 
the accretion rate, and the heating mechanism is indeed ex- 
tremely time dependent. In systems such as Aql X-1 and 4U 
1608-522 the outburst recurrence time is comparable to the 
diffusion time to the surface however. The heat wave result- 
ing from an accretion outburst will therefore be smeared out, 
and the crust temperature will not vary significantly between 
outbursts unless one assumes a high thermal conductivity in 
the crust, so a steady state model using the average accretion 
rate should be able to predict the temperature. Conversely in 
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systems such as Cen X-4 where the recurrence time is signifi- 
cantly longer than the diffusion time to the surface one would 
expect larger temperature variations following an outburst. 

Furthermore in such systems the core - or the surface lay- 
ers of the core in the non-isothermal case - may have time to 
cool significantly between outbursts. In this case we would not 
expect Eqs. (II 2> . (I14> and il6\ to hold, but if (speculatively) 
the tunneling rate late in the cycle is a better predictor of the 
temperature at this point, Eq. il2\ assuming neutrino cooling 
only would actually give a surface temperature around 76 eV 
for Cen X-4 as observed. The same argument for XTE J 1709- 
267 assuming a 6 year recurrence time would lead to surface 
temperatures as low as 83 eV. Here the recurrence time may 
not be quite long enough however, and furthermore the tem- 
perature we cite was measured in May 2003 not very long after 
the March-April 2002 outburst, so the argument may be too 
simple. The temperature we cite for Cen X-4 was measured in 
1997-2000 long after the 1979 outburst so this should provide a 
more clear-cut case. For the temperature of the three remaining 
sources to reach their upper bounds in this way with the cited 
accretion rates, the outburst recurrence times would have to be 
very long (~ 50, 100 and 800 years for SAX J1810.8-2609, 
XTE J2123-058 and SAX J1808.4-3658 respectively). 

This shows how the the time dependence of the heating 
mechanism may potententially help explain the coldest soft X- 
ray transients, even if Qsqm is too high for photon cooling to 
dominate. More work including details on the cooling of the 
core and the transport properties of both crust and core would 
be required though to show that it is possible for the tempera- 
ture to change this rapidly. 

6. Conclusion 

Using a simple scaling argument for the steady state tempera- 
ture of accreting strange stars we have presented a consistency 
check between the constraints on strange quark matter derived 
from superburst ignition conditions, and the temperature of qui- 
escent soft X-ray transients - both assumed to derive from deep 
crustal heating by transfer of matter from crust to core. 

For systems with short recurrence times where the steady 
state temperature may be expected to be a good approxima- 
tion we found reasonable agreement assuming neutrino cooling 
only - although the range of quark matter parameters giving 
superbursters between the isothermal curves in Fig. ^ is then 
excluded by soft X-ray transients. Assuming a combination of 
photon and neutrino cooling for soft X-ray transients, we saw 
in Fig.|3that if Qsqm ^ 1 MeV even the coldest sources may 
be explained while still fulfilling superburst ignition conditions 
at high accretion rates. 

We further showed that the tunneling rate near the equilib- 
rium crust mass depends sensitively on the crust density and 
consequently that essentially all the heating takes place during 
or immediately after an outburst. This points to the need for 
a fully time dependent model to fit systems with recurrence 
times significantly longer than the thermal diffusion time to 
the surface, but if the low level tunneling rate long after the 
outburst is a better predictor for such systems, it would pro- 
vide a fit for Cen X-4 assuming only neutrino cooling as well. 



The upper temperature bounds on SAX J1810.8-2609, XTE 
J2123-058 and SAX J1808.4-3658 can not be reached for the 
cited accretion rates in this way unless the recurrence times are 
very long. This is consistent with the outburst histories of SAX 
J1810.8-2609 and XTE J2123-058 - one known outburst each 
- but does not compare well with the ~ 2 yr recurrence time 
for SAX 11808^4-3658. Similarl y to the sittiation w ith neutron 
stars (iCampana et al]l2002t IVakovlev et al.ll2004l) some form 
of enhanced cooling would thus seem to be required if SAX 
J1808. 4-3658 has accreted for long enough to heat the core to 
its equilibrium temperature and Qsqm is too high for photon 
cooling to dominate. XTE J 1709-267 has an intermediate re- 
currence time, and its temperature was measured shortly after 
an outburst making interpretation difficult. 

If - as seems reasonable - a combination of the two cooling 
mechanisms is required to explain the observations we should 
further note that the constraints imposed by Aql X-1 in Fig.^ 
are not relevant since they were derived assuming only neutrino 
cooling. The requirement that Qsqm ^ 1 MeV can therefore be 
met - even if only just - for the quoted range in Qy and K 
assuming the parameters fit superburst ignition conditions in 
Fig lD Although such fine tuning in the binding energy can not 
be ruled out it is rather conspicuous, but given the simplicity 
of our model and that we have not dealt with the uncertainties 
or the time dependence in detail such quantitative conclusions 
should be treated with caution. We have attempted to test for 
consistency, not to determine the quark matter parameters. This 
would require more detailed modelling, possibly taking the in- 
dividual accretion history of each source into account, which 
may alleviate the need to keep Qsqm small. 

Thus even assuming a high Qsqm the model seems consis- 
tent for the systems where it should be most credible and with 
the exception of SAX J1808. 4-3658 no clear inconsistency can 
be demonstrated for the other sources studied here. However 
SAX J1808. 4-3658 is difficult to explain with neutron stars as 
well, and its precise nature remains to be decided. At this point 
strange stars thus seem to do no worse than neutron stars in ex- 
plaining the quiescent emission of soft X-ray transients regard- 
less of the assumed core parameters, and for a low Qsqm they 
explain even the coldest sources including SAX J1808. 4-3658. 
Because the two cooling mechanisms are important in different 
accretion regimes it finally seems likely that detailed modelling 
along with more accurate temperature measurements at low ac- 
cretion rates would constrain the core parameters further - or 
perhaps demonstrate a clear inconsistency. 

Finally we should not omit a short discussion of the 
interesting qua si persistent transien t s KS 1731-260 and 
MXB 1659-29 ("Wi inands et al J 1200 it iRutledge et al.jEo02bt 
Wijnands et al. 20021 12003L l2004t ICackett et al.ll2006h . These 
are systems with prolonged accretion outbursts - 12.5 and 2.5 
years duration respectively - and a convincing case has been 
built for the suggestion that their crusts are significantly heated 
during these outbursts and brought out of thermal equilibrium 
with the core. The thermal components in their spectrum there- 
fore sho w the cooling of the crust to thermal equilibrium and 
recently fcackett et a l. (2006) using new Chandra and XMM- 
Newton data found that both systems are now close to equilib- 
rium allowing the core temperatures and thermal cooling time 
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scales of the crusts to be determined. Fitting the cooling curves 
with an exponential decay towards a constant temperature pre- 
sumably reflecting the core temperature the constant was found 
to be 71.3 ± 1.6 eV for KS 1731-260 and 52.1 ± 1.6 eV for 
MXB 1659-29 at infinity. Using a canonical mass-radius rela- 
tionship this implies surface temperatures of 93.1 ±2.1 eV and 
68.0 + 2.1 eV respectively, so these sources are both cold but 
not unusually so. Only upper limits to the accretion rates can be 
established with any certainty however due to the very limited 
number of observed outbursts. For KS 1 731-2 60 (one outburst 
in 1988-2001) this is 5.1 x 10"' vr'' dVakov lev et al. 20031) 
using the time-averaged flux of 'Rutledge et aLj ti2002bi') . For 
MXB 1659-29 Wijnands et al. (2003) find bolometric fluxes 
of 5 - 10 X 10"'" erg cm^^ s"' for the 1999-2001 outburst and 
2 - 12 X 10"'° erg cm^^ s"' for the 1976-1979 outburst which 
using the highest of these values in Eq. ( I18> (with / = 1) means 
an upper limit of 1.6 x 10"'" yr"' for a distance of 10 kpc. 
Referring to Fig. |3l both of these limits are very high given 
the surface temperatures, and would require enhanced cooling 
for both strange stars and neutron stars - but then the accre- 
tion rates are hardly certain enough to back up such a claim. 
Perhaps more interesting the e-folding times for the crust tem- 
perature relaxing back to the core temperature were determined 
both sources. Both of these were rather short (246 ± 62 days for 
KS 1731-260 and 501 + 61 days for MXB 1656-29) and by 
comparison wi t h the ne utron star cooling curves computed by 
iRutledge etaP ( l2002bl) for the accretion history of KS 1731- 
260 may be taken as evidence for high thermal conductivity 
in the crust. Strange stars are missing the deep crust between 
10" and 10'^ g cm"^ and as we have seen there is very little 
tunnelUng in quiescence, so although it would go too far to 
show this in detail here, it should be safe to assume that they 
cool faster than neutron stars relaxing the need for high ther- 
mal conductivity in the crust. This new development therefore 
further underlines how rewarding a time dependent model for 
accreting strange stars might be. 
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